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Abstract. The Local Bubble (LB) is an X-ray emitting region extending 
100 pc in radius in the Galactic plane and 300 pc perpendicular to it, and it is 
embedded in a somewhat larger Hi deficient cavity. Its origin and spectral prop- 
erties in UV, EUV and X-rays are still poorly understood. We have performed 
3D high resolution (down to 1.25 pc) hydrodynamic superbubble simulations of 
the LB and Loop I superbubble in a realistic inhomogeneous background ISM, 
disturbed by supernova (SN) explosions at the Galactic rate. We can reproduce 
(i) the size of the bubbles (in contrast to similarity solutions), (ii) the interaction 
shell with Loop I, discovered with ROSAT, (iii) predict the merging of the two 
bubbles in about 3 Myr, when the interaction shell starts to fragment, and, (iv) 
the generation of blobs like the Local Cloud as a consequence of a dynamical 
instability. The O VI column densities are monitored and found to be in excel- 
lent agreement with Copernicus and FUSE absorption line data, showing LB 
column densities < 1.7 x 10 13 cm~ 2 , in contrast to other existing models. 



1. Introduction 

Standard Local Bubble (LB) models fail to reproduce the observed low O VI 
absorption column density (Shelton & Cox 1994; for a recent discussion see 
Breitschwerdt & Cox 2004). Heliospheric in situ measurements are sensitive to 
the boundary conditions imposed by the Local Bubble and the O VI column 
density in absorption is a crucial test for modelling of the local ISM (Cox 2004). 
It seems most plausible that the LB is theresult of 20 successive explosions, 
originating from stars with masses between 11 and 20 Mq from the moving 
subgroup Bl of the Pleiades in the last 14.5 Myrs (Berghofer & Breitschwerdt 
2002). We have simulated the LB evolution by means of a 3D AMR hydrocode 
and calculated the observed O VI column density in absorption (along lines of 
sight crossing Loop I). 

2. Model and Simulations 

We use the 3D model of Avillez (2000), where the ISM is disturbed by supernova 
(SN) explosions at the Galactic rate, and took data cubes of previous runs with 
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a finest adaptive mesh refinement resolution of 1.25 pc (Avillez & Breitschwerdt 
2004). We then picked up a site with enough mass to form the 81 stars, with 
masses, M*, between 7 and 31 M Q , that compose the Sco Cen cluster; 39 massive 
stars with 14 < Af* < 31 Mq have already gone off, generating the Loop I cavity. 

Presently the Sco Cen cluster (here lo- 
cated at (375, 400) pc) has 42 stars to ex- 
plode in the next 13 Myrs). We followed 
the trajectory of the moving subgroup Bl 
of Pleaides, whose SNe in the LB went off 
along a path crossing the solar neighbour- 
hood (Fig. 1). Periodic boundary condi- 
tions are applied along the four vertical 
boundary faces, while outflow boundary 
conditions are imposed at the top {z = 10 
kpc) and bottom (z = —10 kpc) bound- 
aries. The simulation time of this run was 
30 Myr. 
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Figure 1. O VI contour map of a 
3D Local Bubble simulation 14 Myr 
after the first explosion; LB is cen- 
tered at (175, 400) pc and Loop I at 
(375, 400) pc. 



Results 



Morphology: The locally enhanced SN 
rates produce coherent LB and Loop I 
structures (due to ongoing star formation) within a highly disturbed background 
medium. The successive explosions heat and pressurize the LB, which at first 
looks smooth, but develops internal temperature and density structure at later 
stages. After 14 Myr the LB cavity, bounded by an outer shell, which will start 
to fragment due to Rayleigh- Taylor instabilities in ~ 3 Myr from now, fills a 
volume roughly corresponding to the present day size (Fig. 1). 





Figure 2. O VI column density averaged over angles (left panel) indicated 
in Fig. |3 and maximum column density (right panel) as a function of LOS 
path length at 14.1 < t < 15 Myr of Local and Loop I bubbles evolution. 



Oxygen vi: The average and maximum column densities of O VI, i.e., 
(N(0 VI)) and N max (0 VI) are calculated along 91 lines of sight (LOS) ex- 
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tending from the Sun and crossing Loop I from an angle of —45° to +45° 
(s. Fig. 1). Within the LB (i.e., for a LOS length l LOS < 100 pc) (N(0 VI)} 
and N max (0 VI) decrease steeply from 5 x 10 13 to 3 x 10 11 cm~ 2 and from 
1.2 x 10 14 to 1.5 x 10 12 cm -2 , respectively, for 14.1 < t < 15 Myr (Fig. 2), be- 
cause no further SN explosions occur and recombination is taking place. For 
LOS sampling gas from outside the LB (i.e., Ilos > 100 pc) (N(0 VI)) > 
6 x 10 12 and N max (0 VI) > 5 x 10 13 cm~ 2 . The histograms of column den- 
sities obtained in the 91 LOS for t = 14.5 and 14.6 Myr (Fig. 3) show that 
for t = 14.6 Myr all the LOS have column densities smaller than 10 12 ' 9 cm -2 , 
while for t = 14.5 Myr 67% of the lines have column densities smaller than 
10 13 cm~ 2 and in particular 49% of the lines have N(0 VI) < 7.9 x 10 12 cm -2 . 

Noting that in the present model for 
t > 14.5 Myr the O VI column densi- 
ties are smaller than 1.7 x 10 13 cm -2 and 
(N(0 VI)} < 8.5 x 10 12 cm" 2 (blue line in 
both panels of Fig. 2) and that a mean 
column density of 7 x 10 12 cm -2 is in- 
ferred from analysis of FUSE absorption 
line data in the Local ISM (Oegerle et 
al. 2004) we estimate that the present LB 
time is 14.7toJ M Y r - 




LOS length = 100 pc 
Total Number of LOS = 91 

t= 14.5 Myr 
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4. Conclusions 



Figure 3. Histogram of the per- 
centage of LOS with various ranges 
of observed N(0 VI) within the LB 
at t = 14.5 and 14.6 Myr. 



The O VI column density is a sensitive 
tracer of the age of an evolved superbub- 
ble, and can thus give a constraint on the 
timescale of the last explosion, which oc- 
curred in the present simulations 0.61q| 
Myr ago, whereas the nearest explosions to the Sun happened 2.8 ± 0.7 Myr 
ago, in good agreement with the most recent dating inferred from 60 Fe isotope 
analysis in the ferromanganese crust of deep ocean layers (Knie et al. 2004). 
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